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Abstract 

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) 
on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present 
variations of the solar acoustic mode frequencies caused by the solar activity cy- 
cle. High-degree (100 < £ < 900) solar acoustic modes were analyzed using global 
helioseismology analysis techniques over most of solar cycle 23. We followed the 
methodology described in details in Korzennik, Rabello-Soares & Schou (2004) to 
infer unbiased estimates of high-degree mode parameters (see also Rabello-Soares, 
Korzennik & Schou, 2006) . We have removed most of the known instrumental and 
observational effects that affect specifically high-degree modes. We show that the 
high-degree changes are in good agreement with the medium-degree results, except 
for years when the instrument was highly defocused. We analyzed and discuss the 
effect of defocusing on high degree estimation. Our results for high-degree modes 
confirm that the frequency shift scaled by the relative mode inertia is a function of 
frequency and it is independent of degree. 
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1 Introduction 



The correlation between solar acoustic mode frequencies and the magnetic 
activity cycle is well established and has been s ubstantially studied during the 



last and current solar cycles (see for example iDziembowski fc Goodd . l2005l . 
and references therein). The general frequency shift behavior suggests that 
the perturbation by the solar cycle on the acoustic frequencies occurs very 
close to the solar surface. The physical mechanism, however, is still a matter 
of debate, namely whether frequencies are changed directly by the magnetic 



field or indirectly through an associated change in the solar structure (like a 
pressure change). A detailed analysis of the frequency shift characteristics will 
hopefully help understand its physical origin. 

We present here a new analysis of high degree mode frequency shifts over most 
of solar cycle 23, using a global helioseismology technique. High-degree modes 
are best suited for studying the frequency shift correlation with the solar cycle 
since the causing mechanism is believed to be localized near the surface and 
high-degree modes are also confined to a region near the surface. As a result 
the frequency shifts for high degree modes are up to an order of magnitude 
larger than for medium degree modes. 



2 Global Helioseismology Analysis of High-Degree Modes 



In the traditional global helioseismology data analysis methodology, a time 
series of full-disk Doppler solar images is decomposed into spherical harmonic 
coefficients, characterized by its degree {£) and its azimuthal order (m). Each 
coefficient time series is Fourier transformed, and the order of the radial wave- 
function (n) gets separated in the frequency domain. However, a spherical 
harmonic decomposition is not orthonormal over less than the full sphere - 
i.e., the solar surface that can be observed from a single view point- resulting 
in what is referred as spatial leakage. At low and intermediate degrees, these 
leaks are separated in the frequency domain from the target mode and individ- 
ual modes can be identified and fitted. However, at high degrees, the spatial 
leaks lie closer in frequency (due to a smaller mode separation) and, at high 
frequency, they become wider (as the mode lifetimes get smaller), resulting in 
the overlap of the target mode with the spatial leaks that merges individual 
peaks into ridges. The characteristics of the resulting ridge (central frequency, 
amplitude, etc. . . ) do not correspond to those of the underlying target mode. 
This has so far hindered the estimation of unbiased mode parameters at high 
degrees. 

Our methodology to recover the mode characteristics consists in generating 
and fitting a sophisticated model of the underlying modes that contribute to 
the ridge power distribution and dedu ce the offset betw een the ridge proper- 



ties and the underlying target mode (IKorzennikl . Il998l ). To do this we need 



a very good model of the relative amplitude of all the modes that contribute 
to the ridge, i.e. the leakage matrix, which in turn requires a very good 
knowledge of the instrumental properties. The instrumental characteristics 
must therefore be very well under s tood and very precisely measured (see 



Rabello-Soares. Korzennik fc Schoul . I2OOII ). Although, the highly successful 



MDI instrument is responsible for great progresses in our understanding of the 
Sun, there remain some instrumental effects that there are not well character- 
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ized and/or not taken into account in the data analysis (IKorzennik. Rabello-Soares fc Schoul . 
20041 ) ■ 



We have incorporated most of the known instrumental effects important to 
the high-degree analysis in our most recent spatial decomposition, introducing 
one effect at a time to estimate the resulti ng correction on mode parameters 



( iRabello-Soares. Korzennik fc Schoul . l2006l ) . The following known instrumen- 



tal effects were taken into account in this latest spherical harmonic decompo- 
sition, whose results are presented here: (1) the correct instantaneous image 
scale, (2) the radial image distortion, (3) the effect of a tilt of the CCD with 
respect to the optical axis, (4) the effective P angle and (5) a correction to 
the Carrington elements. 

Since we are focusing our attention in this paper on changes of mode param- 
eters with epoch, we are primarily concerned with instrumental effects that 
change over time. Although the MDI instrument has been very stable over 
the years, continuous exposure to solar radiation has increased the instrument 
front window absorption resulting in a continuous small increase of the in- 
strument best focus. Moreover the change of the front window temperature 
due to the satellite orbit around the Sun also adds a small annual variation. 
The instrument has however an adjustable focus with 9 possible positionj^ 
chosen to best suit the science team needs, resulting unfortunately in abrupt 
jumps every time that it is changed. These changes are responsible for the 
largest variations. The best focus is periodically and em pirically determine d 



from intensity images taken with different focus positions (IKuhn et al.l . 120041 ) 



The change in the image scale resulting from a change in focus was shown to 



be an important effect (IKorzennik. Rabello-Soares &: Schoul . |2004| ) and it is 



taken care in the analysis presented here. However, the effect of the amount 
of defocus was not taken into account. A preliminary analysis indicated that 



it ha s a small effect on the frequencies (IKorzennik. Rabello-Soares fc Schou 



20041 ): we have further studied its influence on the ridge frequency. 



3 Description of the Data Set Used 



Full-disk Dopplergrams, acquired by the MDI instrument while operating in 
its 4" resolution mode, were used for this analysis. The MDI instrument is 
operated in this mode, known as the Dynamics Program, every year for two to 
three months, when the required telemetry bandwidth is available. The spher- 
ical harmonic decomposition of these Dopplergrams was computed for every 
tenth I between 100 and 900, and t he resulting time series Fourier transformed 
in shorter segments (as described in lKorzennik. Rabello-Soares fc Schoul . |2004| . 



One focus step corresponds to approximately a third of a wave 
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namely 4096 minutes or 2.8 days long) whose power spectra were averaged to 
produce an averaged power spectruno with a low but adequate frequency 
resolution to fit the ridge while reducing the realization noise. Table 1 lists 
the properties of each epoch included in this study. The peaks in each (£, m) 
spectra were then fitted with asymmetric Lorentzian profiles and the 2i + 1 
frequencies for each (n , i) multiplet expanded in terms of Clebsch -Gordan co- 



efficients up to ae (see iKorzennik. Rabello-Soares &: Schoul . |2004| . for details) 



We also compare our results (hereafter referred to as the high-£ set) to re- 
sults obtained by one of us using the MDI Structure Program. In this mode 
the original Dopplergrams are convolved, on board the MDI instrument, with 
a Gaussian and subsampled on a 200 x 200 grid, thus reducing the teleme- 
try requirements but also limiting the spatial resolution to modes with degree 
£ < 300. The data sets consist of 72-day-long time series that partially overlap 
with the Dynamics periods (see Table 1). The central frequencies and the split- 
ting coefficien ts were determined directly from fit to symmetric Lorentzians 



( ISchod . Il999l ). Although the instrumental effects listed in Section 2 were not 



taken into account here, their effects are m uch smaller because the spatial leaks 



are well separated from the target mode (ISchou &: Rabello-Soaresl . 1200 ll ) 



The coverage of both analyses is illustrate in an i — u diagram in Fig 1. 
Note that there is a small overlap in the coverage of both methods. Indeed, 
by degrading the frequency resolution one can merge resolved modes into 
a smooth ridge and apply the ridge fitting methodology developed for high 
degrees to intermediate degrees. 

The MDI data analyzed here, see Table 1, covers much of solar cycle 23, 
as illustrated in Fig. 2. The average solar UV spectral irradiance (given by 
the NOAA Mg II core-to-wing ratio) relative to the 30-day running mean 
maximum of solar cycle 23 (reached in January 2002) for each epoch is listed 
in Table 1. That table also lists the average defocus of the MDI images during 
each Dynamics epoch, in units of MDI focus steps. Early on (1996 to 1998) 
the instrument was set to a large defocus (up to 2.5 focus steps), while later 
on the instrument was set to be nearly on focus {i.e., a defocus of 0.5 focus 
steps in average). 



4 Analysis of the Frequency Shift 



Figure 3 shows mode frequency differences between 2002 and 2005, namely be- 
tween high and moderate activity epochs respectively, resulting from medium 

^ Henceforth the number of averaged spectra varies with epoch, from 12 segments 
in 2003, to 32 in 1997. 
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and high degree determinations. The frequency shifts for each order, n, scale 
as a power law with the quantity v/L (where L = I + 0.5, a good proxy 
of the mode inner turning point), and can be described in a first approxi- 
mation by Sfn/ oc 1^], wh ere In/ is the mode inertia (as first noticed by 
Libbrecht fc Woodardl .' 1l99oh . 



For the observational periods where the MDI instrument was set to a large 
defocus (1996 to 1998), our measured frequency shifts for high-£ do not agree 
with results from medium-^. The defocus has a substantial effect on the ridge 
central frequencies, contrary to our initial expectations. The resulting instru- 
mental point spread function (PSF) is affected and the effect of a varying PSF 
has not yet been incorporated in our analysis, therefore we will not include 
these epochs in our analysis of changes of mode parameters with the solar 
cycle. 

Figure 4 shows the correlation between the frequency shift and the solar ac- 
tivity index variation (given by the solar UV irradiance) for four modes with 
frequency around 3 mHz at high degree (^ = 200, 400, 600 & 800) for the 
1999 to 2004 epochs with respect to the 2005 epoch. We have fitted, using 
a weighted least-squares minimization, the frequency shift to the solar index 
variation, assuming a linear relationship with a zero intercept (solid line). 
This fit was performed for each (n, ^) mode of the medium and high-^ sets 
separately, resulting in a total of 1468 and 423 regression fits respectively. 

We also computed correlations between the frequency shift and the Magnetic 
Plage Strength Index (Mt. Wilson Observatory), the solar radio 10.7 cm flux 
(National Research Council of Canada) and the sunspot number (SIDC, RWC 
Belgium) using the high-£ set. The magnetic index gives similar Pearson cor- 
relation coefficients as the UV flux while the radio flux coefficients are slightly 
smaller (by a few percent). The sunspot number correlation coefficients are 
also similar to those for the UV flux but only for frequencies around 3 mHz, 
their correlation coefficients increase as the frequency decreases, by some 20%. 

The frequency shift between solar maximum and minimum can be estimated 
by multiplying the slope of the linear fit (solid line in Fig. 4) by the corre- 
sponding solar index variation. We used the 30-day running mean values to 
estimate the minimum and maximum, reached in March 1996 and January 
2002 respectively. We rejected modes whose Pearson correlation coefficient is 
smaller than 0.8, or whose slope uncertainty is larger than 20% and only fitted 
modes for which we have estimates for at least five of the six epochs. This es- 
timated minimum-to- maximum frequency shift, 5z/^, incorporates information 
from all the data sets (1999 to 2005) in a consistent way, decreasing ipso facto 
the uncertainties. 

The frequency shift increases sharply with frequency and degree. Although 
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several authors multiply the frequency shift by the mode inerti a to remove its 



dep endence, we found, as others have at lower and medium-^ (j Chaplin et al. 



20011 ). that the mode inertia normalized by the inertia of a radial mode of the 
same frequency gives a better agreement between the medium and high-£ set 
shifts, as shown in Fig. 5. The estimated minimum-to-maximum frequency 
shifts can be described according to a simple power law, i.e.: 

= %^ (1) 



where Qn,i is the relative mode inertia, i.e.: 

Qn,e = j-^ (2) 



and 
S 



see 



is the mode inertia calculated fr om Christensen-Dalsgaard's model 



Christensen-Dalsgaard et al.l . Il996l ). The frequency dependence of the 



frequency shift scaled by Q„ £ shows the relative importan ce of the contribution 



from the surface layers in relation to the the solar interior ( IChristensen-Dalsgaard et al. 
1989h . 



We fitted Eq. 1 to estimate 7, using a weighted least-squares minimization, 
as illustrated in Fig. 5. For p modes, we found very similar results for the 
medium- and high-£ sets: 7p = 3.64±0.02 (dashed) and 7p = 3.62 ±0.04 (solid) 
respectively. Note that in the lower panel there is a step in the frequency shift 
around 2.3 mHz and only modes with u > 2.5 mHz were included in the fitting 
of the medium-^ set mentioned above. Despite the step, modes with z/ < 2.23 
mHz show a similar value of 7 as for i> > 2.5 mHz: 7^ = 3.5 ± 0.3 (dotted and 
dashed line). It is interesting to point out that the height of the upper turning 
point decreases very fast with mode frequency for z/ < 2.3 mHz (see Fig. 2 in 



Chaplin et al.l . I2OOII ). The / modes for the high-£ set are strongly correlated 



with the solar cycle (see lower right panel in Fig. 4) and exhibit a behavior 
similar to the p modes with a positive 7/ but with a 7/ = 1.6 ± 0.2 (dotted 
line in the upper panel of Fig. 5), less than half of 7^, which could be an 
indication that diffe rent physical effects are respon sible for the solar- induced 
frequency shifts (see iDziembowski fc Goodd . 120051 ). There are no / modes in 
the medium-^ set with a correlation higher than 0.8 for comparison. 



Chaplin et al.l (I2OOII ) scaled the frequency shift by the mode inertia, i.e.: 



Ca 



i-nl 



(3) 



and estimated a for 4 < £ < 150 using GONG data. They found a ~ for 
1.6 < z/ < 2.5 mHz and a = 1.92 ± 0.03 for 2.5 < u < 3.9 mHz. To compare 
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our results, we multiplied our fit, u"^, by In,o and fitted the function Ca z^"- 
Our results using MDI medium and high-degree modes agree well in the high- 
frequency range, i.e.: a = 1.97 ±0.02. However, MDI medium-^ modes exhibit 
a negative slope a = —3.05 ± 0.01 for z/ < 2.5 mHz. Note that log((5z/^£ In/) 
display a quadratic dependence on log(z/„^£), with an inflection point at 2.59 
mHz as shown in Fig. 6. This minimum in the frequency shift scaled by th e 



mode inertia In/ is also seen in Figures 4 and 5 in lHowe. Komm fc Hill! (120021 ) 



The solar-cycle correlation coefficients decrease sharply for modes with n > 6 
and frequencies larger than ~ 4.8 mHz, some of them becoming even negative. 
However, these "anti-correlation" coefficients are not large enough to draw any 
conclusion other than these modes present no significant correlation with the 
solar cycle. 



5 Influence of Image Defocus on the Frequency Shifts 



As mentioned before, the frequency shifts for 1996 to 1998 are biased due to the 
change in the amount of defocus. To estimate the resulting bias, we computed 
the correlation versus the solar activity index and estimated the minimum- 
to-maximum frequency shift, Su'^, as described in the previous section, but 
using only these three epochs. These biased frequency shifts remain strongly 
correlated with the solar cycle, but the result is affected by the defocus, as 
illustrated in Fig. 7. The differences are substantial for n < 4 and for i > 600. 
Note that the change in the amount of defocus with respect to 2005 averages 
to 0.16 (in focus steps) for the later years (1999 to 2004), while that same 
quantity averages to 1.77 for the early years (1996 to 1998), namely an order 
of magnitude larger. The remaining frequency bias due to change in defocus 
could be as large as the frequency shift due to solar activity. 

We can not discard the fact that the frequency shifts in Fig. 5 using 1999 to 
2005 may have some significant, although very small, bias due to variations 
in the defocus between the data sets {i.e. 1999 to 2004 with respect to 2005). 
The instrumental PSF must be taken into account in the analysis at least for 
data sets taken when the instrument was out of focus. 



6 Conclusions 



We analyzed the variation of the solar acoustic frequencies during solar cycle 23 
using medium and high degree modes. Our high-£ mode analysis incorporates 
most observational and instrumental effects that affect specifically high-£ in 
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the spherical harmonic decomposition to produced unbiased mode estimates 
from ridge fitting. 

The frequency shift scales very well with the mode inertia normalized by 
the inertia of a radial mode of the same frequency and follows Eq. 1, where 
7 = 3.63 ± 0.02 for p modes using medium-^ (MDI Structure Program) and 
high-£ (MDI Dynamics Program) data and 7 = 1.6 ± 0.2 for the / modes, 
which exhibit a different behavior. 

The frequency variation does not appear correlated with solar activity for 
frequencies larger than ~ 4.8 mHz. There is also an abrupt step in the scaled 
frequency shift for frequencies smaller than ~ 2.3 mHz. 

We also analyzed the effect of defocusing on the high-^ frequency determina- 
tion. This is particularly important for the first three years of MDI observa- 
tions when the instrument was set to a substantial defocus. This has a larger 
impact on modes with n < 4 and i > 600. Further progress will be achieved 
by including a time varying PSF in the high-£ global helioseismology analysis. 
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Table 1 

Observational periods, relative average solar active index and average defocus. 



Year 


Starting Date 

Dynamics 
[Month/Day] 


Duration 
Dynamics 
[days] 


Solar Index 
rel. to max. 

[in %] 


Average 
Defocus 


Starting Date 

Structure 
[Month/Day] 


1996 


05/23 


63 


2±3 


2.07 ±0.07 


05/01 


1997 


04/13 


93 


5±3 


2.53 ±0.07 


04/26 


1998 


01/09 


92 


20±7 


1.75 ±0.05 


02/08 


1999 


03/13 


77 


40 ±13 


0.1±0.4 


02/03 


2000 


05/27 


45 


69 ±11 


0.45 ±0.02 


04/10 


2001 


02/28 


90 


61 ±14 


0.57 ±0.04 


01/23 


2002 


02/23 


72 


80 ±8 


0.77±0.05 


03/31 


2003 


10/18 


38 


51 ±17 


0.92 ±0.01 


10/28 


2004 


07/04 


65 


36 ±11 


0.2±0.1 


08/11 


2005 


06/25 


67 


30±9 


0.34 ±0.02 


05/26 



5000 




200 400 600 800 

Spherical Harmonic degree I 



Fig. 1. Coverage, in an £ — diagram, of the medium-^ (crosses) and high-^ (circles) 
set for 1996. The coverage is very similar for all epochs. Modes with I < 200 and 
v < 4 mHz may not completely overlap with their spatial leaks and were excluded 
from the present high-^ analysis. 
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Fig. 2. Solar UV spectral irradiance (given by the NOAA Mg II core-to-wing ratio) 
30-day running mean relative to its maximum during solar cycle 23 (line) and its 
average for each epoch. 
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Fig. 3. Differences between the mode frequencies as observed in 2002, near solar 
maximum, and in 2005, during moderate activity. Only modes with n = to n = 8 
are plotted (evens in black, odds in grey). The medium-^ set is indicated by stars 
(lower part of the plot) and high-£ set by circles (upper part). 
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Fig. 4. Examples of correlation between the frequency shift and the solar cycle for 
six epochs (1999 to 2004) with respect to 2005 for four high-£ modes with different 
degree and order but similar frequency (~ 3 mHz). The error bars are the standard 
deviation of the fitted frequencies. The straight line is a weighted least-squares fit 
to the data, keeping the intercept to zero. 
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2000 3000 4000 5000 




Fig. 5. Minimum-to-maximum frequency shift multiplied by the normalized mode 
inertia as a function of frequency. The upper panel shows only the high-^ set and, in 
the lower panel, the high-£ p modes are indicated in black (circles) and the medium-^ 
p modes in gray (crosses). 
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Fig. 6. Minimum-to-maximum frequency shift multiplied by the mode inertia as a 
function of frequency. The medium and high-^ p modes are represented by crosses 
and circles respectively. The dashed line corresponds to the dashed line in the 
lo wer panel of Fig. 5. The solid line is given by Eq. 3 using the values estimated 
bv lChaphn et al.l (|200lh : a = for 1.6 < i/ < 2.5 mHz and a = 1.92 for 2.5 < u < 3.9 
mHz; the fitted values for Ca were not given by the authors hence ad hoc values 
were used. 



71=0 71=1 n=2 




Fig. 7. Minimum-to-maximum frequency shift estimated using 1999 to 2004 data 
sets with respect to 2005 (black circles - same as in Fig. 5) and using 1996 to 1998 
data sets with respect to 2005 (solid line). Note that, for the later, there is a larger 
scatter and the standard deviation is also larger (not plotted for clarity) due to the 
small number of available epochs. Modes whose Pearson correlation coefficient is 
smaller than 0.8 or whose slope uncertainty is larger than 20% were excluded from 
the plot. 
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